Probing the power of an electronic Maxwell Demon 
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We suggest that a single-electron transistor continuously monitored by a quantum point con- 
tact may function as Maxwell's demon when closed-loop feedback operations are applied as time- 
dependent modifications of the tunneling rates across its junctions. The device may induce a current 
across the single-electron transistor even when no bias voltage or thermal gradient is applied. For 
different feedback schemes, we derive effective master equations and compare the induced feedback 
current and its fluctuations as well as the generated power. Provided that tunneling rates can 
be modified without changing the transistor level, the device may be implemented with current 
technology. 



Maxwell's demon - a hypothetical intelligence in a box 
\ capable of sorting hot (fast) and cold (slow) atoms to 
i left and right sub-cavities simply by matchingly insert- 
i | ing and removing an impenetrable wall - was initially 
^ . suggested by J. C. Maxwell to highlight that thermody- 
namics is a macroscopic theory. Under ideal conditions, 
inserting and removing the wall would not require work 
and the apparent contradiction with the second law - 
the entropy in the box would be reduced and after the 
sorting process, work could be extracted from the tem- 
perature difference between the sub-cavities - has been 
a source of much debate ever since 1 . It is now generally 
believed that this paradox is overcome by the Landaucr 
principle 2 : Deleting the data required for the process- 
ing in the demons mind would at least generate entropy 
S = ks hi 2 or dissipate heat of at least Q = ksT In 2 per 
bit of information, yielding a net production of entropy. 

The demon performs a measurement on the system (is 
the atom slow or fast) and conditioned on the measure- 
ment result it performs an action (opening or closing the 
wall), which is formally nothing but a closed-loop feed- 
back control scheme. In our proposal, wc replace the two 
sub-cavities of the box by two conductors that we assume 
to be in separate thermal equilibria. The conductors are 
coupled by a single resonant level, which does not require 
charging and spin 3 effects. Conditioned on its occupation 
dynamics - provided by a nearby quantum point contact 
(QPC) - the tunneling rates to the conductors are mod- 
ified in time. The feedback schemes we consider here are 
illustrated in Fig. 1. We will address the power produc- 
tion, current and statistical properties of the device. 



I. METHOD 

A single-electron transistor (SET) coupled to two leads 
in thermal equilibria may in sequential tunneling (weak 
coupling) approximation be well described by a Marko- 
vian Lindblad-type master equation p(t) = Cp(t), where 
- when written as a vector - p(t) = (po{t), Pi (t)) de- 
scribes the SET part with occupation probability pi(t). 
The Liouvillian C can be derived from microscopic cal- 
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FIG. 1: (Color Online) Sketch of the setup: An SET with 
two attached leads (sketched are their Fermi functions /l 
and fn) is capacitively coupled to a nearby QPC. Its current 
Iqpc yields information on the instantaneous occupation of 
the SET. One may either apply Liouvillians jCe and Cf with 
different tunneling rates conditioned on the present state of 
the SET (scheme I, compare different tunneling barriers in 
SET sketch) or control operations ki and no conditioned on 
the change of the SET state (schemes Ila/b, not shown in 
SET sketch). During the control operations, electrons may 
also tunnel into or out of the SET (dashed shorter spikes), 
which may not (Ha) or may (lib) recursively trigger further 
control operations. 



culations or simply using Fermis golden rule and reads 4 
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where T Q denotes the tunneling rate between SET and 



lead a and f a 
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the Fermi func- 



tion of lead a with corresponding inverse temperature 
j3 a and chemical potential p a at the SET level e (as- 
sumed to be constant throughout). Since the total 
number of electrons is conserved in the tunneling pro- 
cesses, it is possible to uniquely identify matrix ele- 
ments of the Liouvillian with electronic jump processes 
into and out of the left and right leads, which en- 
ables one to convert Eq. (1) into an infinite set of 
coupled equations for (ul, n#)-resolved density matrices 
p{n L ,n R ) ^ _ These are conditioned on the number of elec- 
trons that have tunneled after time t into or out of the 
left (til) and right (tir) leads. Fourier-transformation, 



p(Xl,Xr) = Hn L ,n R P 



(n L ,n R ) n\ e iXLn L p ix 



e ,XanR , yields the 
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compact equation p(xL,XR,t) = £{XL,XR)p(XL,XR,t) 
with 

£(XL,XR) - a ^ R} Ta [+f a e-^ _(i_ /a ) J 

= £ r Q j- Q ( XQ ). (2; 

a£{L,R.} 

In absence of feedback, the above equation yields the 
complete statistics of electrons (current, noise, skewness 
etc.) that have tunneled across the SET 5 . At small 
SET bias, the QPC does not resolve to which of the at- 
tached leads a tunneling process has happened. There- 
fore, simplest feedback operations can only be condi- 
tioned on whether the SET is empty or filled (scheme 
I) or whether an electron has jumped into or out of the 
SET (schemes Ila/b). Such control operations may easily 
be included in the numerical solution of Eq. (2) via an as- 
sociated stochastic Schrddinger equation 6 . For analytic 
results however, it is more favorable to derive an effective 
Liouvillian under feedback (compare appendices A, C, 
and D), which enables one to systematically study the 
effects of feedback on the statistics. 

In scheme I, wc only apply two different Liouvil- 
lians - conditioned on whether the SET is empty 

(£e{xl,Xr) = Y, a e{L,R} T a eSaj: a (Xa)) or filled 

{C-f(xl,Xr) = Y, a e{L,R.} r " e a F a {x a )), where the di- 

e If 

mcnsionless feedback parameters S a ' e M encode the 

~E IF 

modification of tunneling rates and 5 a = recovers 
the situation without feedback. The effective feedback 
generator reads in this case 

£fb(X£, Xr) = £e(xl,Xr) (q Q^j 

+£f(xl,xr) ( 1 ) • ( 3 ) 

An extremal form of this feedback scheme would be 
to cut the left junction as soon as the SET is filled 
(5?[ — > —00) and to cut the right junction when it is 
empty (5^ —> —00), which automatically implies unidi- 
rectional transport - independent of potential or temper- 
ature gradients. This effectively implements a feedback 
ratchet with two teeth - compare SET tunneling barri- 
ers in Fig. 1. However, in the idealized classical limit 
where the electrons are localized either on the left lead, 
the SET, or the right lead, this ratchet does not directly 
perform work as its potential does not change where the 
electron is localized. The resulting effective Liouvillian 
cannot generally be written in the form of Eq. (1) us- 
ing modified tunneling rates r Q > and Fermi functions 

L e (0,1). 

In scheme Ha, we instantaneously modify the SET tun- 
neling rates by (5-kicks (compare appendix B) immedi- 
ately after an electron has jumped in (e Kl ) or out (e K °) 



of the SET, such that 

ki{xl,Xr) = fiiFaiXa) , 

a£{L.R\ 

ko(xl,Xr) = Y S °F<x{Xa), ( 4 ) 

a£{L.R\ 

where S 1 ^^ > are dimensionless parameters roughly 
given by the product of pulse width and height of time- 
dependent SET tunneling rates (compare appendix B). 
The counting-field dependence arises from the simple fact 
that during the control operation electrons may tun- 
nel through the junctions (in scheme Ha, these tunnel- 
ing events do not trigger further control operations): In 
fact, for infinitely strong feedback at one junction and 
zero feedback at the other junction (e.g., <5£ — > 00 and 
<5jj — » 0) one obtains immediate equilibration of the 
SET with the lead to which it is tunnel-coupled (e.g., 
e «/(0,o)p _^ (1 — /l, /l) for all p). These tunneling events 
have to be taken into account when the complete statis- 
tics is required. Using from Eq. (2) the decomposition 

T a ( Xa ) = J^ + T+e+^ + J^e-*** , (5) 

where (J~~) is responsible for jumps into (out of) 
contact a and conserves the system occupation, the 
full effective feedback Liouvillian reads 

C^ a {XL,XR) = E ^ 
a£{L,R} 

a£{L,R.} 

+e *i(xL,xn) T a T-e-^ ,(6) 

a£{L,R.} 

such that jumps in and out of the system are immedi- 
ately followed by the respective control operation. In con- 
trast to feedback schemes modifying the system Hamil- 
tonian 7,8 it is no longer possible to simply shift counting 
fields by control operators. 

For an infinitely fast QPC sampling rate, it would 
also be possible to recursively trigger further control op- 
erations when tunneling events take place during con- 
trol. In scheme lib we restrict ourselves to the case 
that during the control operations only one junction ad- 
mits tunneling at a time, i.e., 5 ! L = 5^ = 0. Then, we 
may also derive an effective feedback master equation 
for a recursive feedback scheme, where electrons tunnel- 
ing during a control operation would induce further con- 
trol operations - possibly ad infinitum (at infinite bias 
and infinite feedback strength). In this case, the con- 
trol operations only depend on a single counting field 
and we may use the decompositions (see appendix D) 
e *o(xL,XR) = V°(xl) = T>n + T ? 8 e+iXL + Vfe~ lXL 
and e K HxL,x«) = V 1 {xr) = V T N + V\er 1 ^ + V I e +l ^ R 
together with the evident relations Vq(J-^ + J-^) = 
V\[Tl +J r i)= VSK = v \ v ? = ( these effectively 
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imply that during the control operations, the transport 
is unidirectional) to sum up the infinitely many terms 
as a von Neumann operator series. Eventually (see ap- 
pendix D), this results in the effective feedback Liouvil- 
lian 

C^ b ( X L,XR) = 




N''0) 



1 — pO'p* ) e i ( XR ~ XL ^ 



(7) 



Eqns. (3), (6), (7) yield the complete statistics for the 
current through the SET under the different feedback 
schemes: The generating function for its cumulants is in 
the long term given by the dominant (with A(0, 0) = 0) 
eigenvalue \(xl,Xr) 01 the effective Liouvillians. We 
have checked numerically (see appendixes A, C, and D) 
that the analytic results from the effective master equa- 
tion coincide for all schemes with an ensemble-average 
of the stochastic Schrodingcr equation with feedback ex- 
plicitly included. In the idealized classical limit, none 
of the schemes performs work on the system. However, 
quantum-mechanically the time-dependent modification 
of the tunneling rates changes the energy spectrum of the 
total Hamiltonian and thereby performs work 9 . 



II. CURRENT AND FLUCTUATIONS 

We summarize the behavior of the current for finite 
feedback strengths at reverse infinite bias (Jl = 0, 
fa = 1), zero bias (/_ = fu = /), and infinite bias 
(/_ = 1, Jr = 0) for the different feedback schemes in 
table I. For finite temperatures (where < / < 1) feed- 
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TABLE I: Values of the current (in units of T L = F a = T) 
under finite but symmetric feedback. The feedback parame- 
ters have been chosen as 8^ = 5^ = — <5fb and 6r = 5® = <5fb 
for scheme I and S' R = $2 = Sa> as well as 8 ! L = 5§ = for 
schemes Ila and lib, respectively. 



device acts as a demon shuffling electrons from one bath 
to another (here from left to right) using only information 
on SET occupancy. This has to be contrasted with feed- 
back ratchets 10 where the time-dependent ratchet poten- 
tial performs work. The apparent divergence of the cur- 
rent in feedback scheme I for infinite bias and feedback 
strength is natural as the tunneling rates diverge likewise. 
For feedback scheme lib however, we observe a genuine 
feedback catastrophe for large bias and large feedback 
strength leading to a divergent current: Control oper- 
ations mutually trigger further control operations with 
vanishing halting probability, which leads to avalanche- 
like transport. Therefore, in contrast to scheme I, the 
Fano factor F = S/\I\ given by the ratio of noise and 
current - as summarized in Table II - also diverges in 
this region when we let the feedback strength go to in- 
finity. 
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1/2 

4-3e~' 5 fb 


(l-2/)^[cosh(i fb )-l]+sinh(«5 fb ) 
2e 25 fb-e a fb+l 
e 24 fb-l 


4_2 e - a fb 
e Slh - 1/2 



TABLE II: Values of the Fano factor (assuming <5fb > 0) for 
zero and infinite bias with the same parameters as in table I, 
the zero-bias Fano factor for scheme lib has been evaluated for 
/ = 1/2 only for brevity. The reverse infinite bias Fano factor 
(not shown) is 1/2 for all schemes independent of feedback. 
Divergence of the infinite bias Fano factor for scheme lib at 
infinite feedback strength 5fb — ► oo demonstrates the feedback 
catastrophe. Also, at zero feedback the Fano factor should 
diverge at zero bias as this is the point where conventionally 
the current would vanish. 

Naturally, there exists a parameter regime where the 
device transports electrons against an existing electri- 
cal or thermal gradient, where e.g., electrons are trans- 
ported from left to right even though /_ < fn, compare 
also Fig. 2. For scheme I, this effectively implements a 
Parrondo game (ratchet): Playing two losing strategies 
(tunneling with the bias) in an alternating manner may 
yield a winning strategy (tunneling against the bias) 11 . 
For our model, feedback is necessary to achieve such 
a current inversion: Without feedback, the long-term 
cumulant-generating function for the current obeys the 

analytic relation A(0, -\) = A (o, +x - ihi j^^rf^ 
which when both leads are at the same temperature even 
tually leads to the fluctuation theorem 12-14 1im Pn ^> 



back may induce a current even at zero bias, such that the 



e +np(n L -ii R ) _ xhis implies that for a constant Liouvil- 
lian (1), the current / = lim t _j.ee d/dt nP n (t) always 
flows from reservoirs with large chemical potential to- 
wards the reservoir with the smaller chemical potential. 
Also with open-loop feedback, where control operations 
are applied unconditionally (e.g., in a random or periodic 
manner) , we find that the current will never flow against 
the electro-thermal gradient - as long as no work is per- 
formed on the SET itself (i.e., the SET level e remains 
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unchanged), see appendix E. 

For scheme I we find that the Johnson-Nyquist rela- 
tion normally relating noise with differential conductance 
at zero bias voltage (at similar temperatures left and 
right) is now shifted to the equilibrium voltage V* 

Hf - i ~ 

I(V*) = 0, i.e., 

§ Iv I v— v ' ^ e obtain a similarly simple modification of 
the fluctuation theorem for the electron counting statis- 
tics under feedback: The cumulant-generating function 
for the current obeys 



S£ + s F 



L i I (3 at which the current vanishes 
more explicitly we have S\ v=v , = 



\ I (-X) = X I +X-*m 



(1 - h)f. 
h(l-fR 



(8) 



which for both leads at the same temperature yields 
for the counting statistics lim rr^m — e n ^- L ~ At - R_v \ 

t— too 

This formula directly demonstrates that the current may 
be directed against a potential gradient when the feed- 
back parameters and hence V* are adjusted accordingly. 

For feedback schemes Ila and lib, the equilibrium volt- 
age V* can be obtained numerically and depends on the 
baseline tunneling rates F a . 



III. MAXIMUM POWER AT MAXIMUM 
FEEDBACK 
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FIG. 2: (Color Online) Current-voltage characteristics under 
feedback strength 5f, = 8r = — <5j| = — <5£ = 1 in scheme 
I (thin solid curve, red), maximum feedback in schemes Ila 
and lib (dashed and dotted, respectively) and no feedback 
(bold black) for symmetric tunneling rates Yl — Y — Yr and 
fie = 2. The current may point in the other direction than the 
voltage leading to a positive power generated by the device. 
The inset (for fie/C = 1) shows the mean-field evolution of 
the voltage from an initial value (empty circles) to the stable 
fixed points (filled circle) for symmetric capacitances of the 
two conductors. 



and S ! L 



For schemes Ila and lib we may take 8^ = 8° — > oo 
8^ = to obtain the maximum effect, but 
to bound the Liouvillian for scheme I we constrain our- 
selves to finite feedback strengths 8f — 8^ — +<5fb and 
= = — <5fb- The power P = —IV generated 
by the demon can - when both leads arc at the same 
temperature /3 be evaluated as P/3 = —I/3((j,l — (J-r) = 



-I(f L ,f R )\n 



fUl-fn) 



The right hand side of this ex- 



pression can for given and Y R be numerically maxi- 
mized with respect to the Fermi functions /l/r £ (0, 1) 
yielding the maximum power generated by the device, see 
table III. The maximum power generated with scheme I 



scheme 


Y L = Y R = Y 


Tmax Tniin 


Ila 


0.0841 rfc s T 


0.1219 r max fc B T 


lib 


0.1892 TfcsT 


0.1589 r max fc s T 



TABLE III: Maximum power generated with the device for 
schemes Ila and lib at maximum unidirectional feedback for 
symmetric and highly asymmetric tunneling rates. 

is unbounded (as the Liouvillian) and asymptotically ap- 
proaches 0.2785 Te Slb kBT for symmetric tunneling rates. 
Even under idealized conditions, the associated work per- 
formed by the demon has to be contrasted with the heat 
dissipated when the QPC trajectory's data-points are 
deleted (Landauer principle). To perform the continu- 
ous monitoring of the SET it is important that the QPC 
sampling rate At" 1 is greater than the maximum tun- 
neling rate At < T"^ (or At < r max e -5 for scheme I). 



This implies that the maximum work per current mea- 
surement W = P max Ar is always smaller than the dissi- 
pated Landauer heat Q > ksT\n2 by the demon, such 
that the results are compatible with the second law. 



IV. FEEDBACK CHARGING EFFECTS 

Small leads of finite capacitances will usually be driven 
out of equilibrium due to transport. Additional larger 
reservoirs at thermal equilibrium may however immedi- 
ately re-enforce equilibrium in the leads solely by scat- 
tering interactions (without electron tunneling). We may 
phenomcnologically include this effect by making the 
chemical potentials dependent on the number of tun- 
neled particles: Starting at equilibrium, the difference 
V in chemical potentials between the two leads is clas- 
sically simply proportional to the number of tunneled 
particles V = en/C, where C denotes the total capac- 
itance (similar relations hold for left and right chem- 
ical potentials in case of asymmetric capacitances 15 ). 
We may numerically solve the resulting mean-field non- 
linear (compare appendix F) differential equation V — 
e/C (h(t)) = eI(V)/C to obtain the dynamical evolu- 
tion V(t), compare the inset in Fig. 2. After the equilib- 
rium voltage (/(V*) = 0) has been approached, feed- 
back may be stopped and the current will reverse its 
direction (black trajectory). With the QPC already in 
place, it appears reasonable to use it as a detector to 
clearly discriminate the resulting initial fluctuations from 
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equilibrium ones. This scheme works as an accumula- 
tor undergoing charging (feedback) and discharging (no- 
feedback) cycles, where the total stored energy is given 
by W = C{V*) 2 /(2e). 

V. EXPERIMENTAL IMPLEMENTATION 

With gate voltages, the height of the tunneling bar- 
riers may be adjusted such that the timescalc on which 
electrons tunnel through the QPC can be tuned from 
ms to seconds and even hours. For periodically varying 
gate voltages, the frequency corresponds with 100 MHz 
to switching times five orders of magnitude smaller in re- 
cent electron pumping experiments 17,18 . The bandwidth 
of typical experimental QPC detector devices has been 
reported in the range of 40 kHz, with sufficiently larger 
current sampling frequencies of 100 kHz 19 . The experi- 
mental challenge therefore clearly lies in the necessity of 
strongly modifying the tunneling rates without perform- 
ing work on the system (changing the SET level). With 
gate electrodes of sizes below 100 nm 18 it should be pos- 
sible to keep the energy level of the SET (size about 300 
nm 16 ) approximately constant. 

VI. SUMMARY 

To conclude, we have compared several closed-loop 
feedback schemes implementing Maxwell's demon by 



means of an effective feedback master equation. Scheme 
I used a piecewise constant Liouvillian conditioned on 
the time-dependent SET occupation, whereas schemes 
Ha and lib were conditioned on its change and used S- 
kicks in the tunneling rates. All schemes are capable of 
generating a current against a moderate bias - which is 
for constant SET level not possible for open-loop schemes 
- and may for finite-size leads be used to charge a feed- 
back battery. With Landauers principle, the second law 
is of course respected by the device. Schemes I and Ha 
(with necessarily small feedback strength) should be im- 
plementable with present-day technology, whereas the 
feedback recursion depth appears to be currently lim- 
ited by the QPC sampling frequency. Schemes with fi- 
nite recursion depth are however also treatable with the 
methods in this article. 
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Appendix A: Justification of scheme I 



Assuming that the dot is empty at time t, i.e., 
p(t) = (1,0), its no-measurement time evolution un- 
der feedback scheme I will be governed by p(t + At) = 
e c E (xL,XR)^t p(t), whereas for an initially filled dot 
p{t) = (0,1), we will have the evolution p(t + At) = 
e c F (xL,XR)^t p{t). Using projection superoperators on the 
empty and filled dot states, respectively, both cases can 



6 



be incorporated into a single equation 
s £s(xz,,XH)At (10 



p(t + At) = e L 







Pit) 



for the finite-difference [pit + At) — pit)]/ At yields the 
effective feedback Liouvillian Eq. (3) when At — > 0. The 
switching between the different propagators is assumed 
to be instantaneous, such that during the switching time, 
no particles may tunnel. 



+e 




f 



pit). (Al) 



Expanding the propagators for small At, using that the 
projection operators add up to the identity, and solving 



J 



Therefore, the original counting fields may be kept, and effectively, the Liouvillian (3) has the first column from 
Ce and the second column from Cp 



£kiXL,Xn) = 



-T L e s Lf L -T R e 5 *f R 



^ e-^f L + T R e s *e- i **f R 

which explicitly breaks detailed balance. 



-T L e s ^e +i ^(l - / i )+r fl e^e+ i ^(l - f R ) 
-T L e s ^0--f L )-T R e s n(l-f R ) 



(A2) 



r 



A single trajectory for this feedback scheme may 
be generated as follows: Assuming that e.g., the dot 
is initially filled, the probability that during a small 
timcstep At the electron jumps out to the left is given 
by PxeftCAt) = AtT L e s t{l - f L ), whereas the prob- 
ability to jump out to the right reads -P r f ght (At) = 

AtT R e R (l — f R ). The probabilities for electrons jump- 
ing from left or right lead into an initially empty dot are 
obtained similarly and read P^ it (At) = AtT^e L Jl and 
Pj? ht (Ai) = AtT R e SR f R . Using a random number gen- 
erator, one may with sufficiently small timesteps (such 
that the jump probabilities are significantly smaller than 
one) generate single trajectories for the dot occupation 
nsETj the number of tunneled particles to the left lead 
«i, and the number of particles tunneled to the right 
lead n R . The ensemble average of many such trajecto- 
ries may now be compared with the analytic solution of 
the effective feedback master equation (3), see Fig. 3. It 
follows that for sufficiently many trajectories the same 
average observables will be obtained as with the effective 
feedback master equation. 



Liouvillian. 

Let us therefore assume a constant baseline Liouvil- 
lian superimposed with an additional control Liouvillian 
Lit) = Cq + £ c 9(i — i c )0(i c + t c — t), where the latter 
is turned on at time t c and lasts for timespan r c (O(x) 
denotes the Heaviside step function) . The solution for all 
times reads in this case 



Pit) 



e Cot p 



e (C +C c ){t-t a ) e C t c 



Po 



e Co(t-t c -r c ) e (Ca+C a )r c e £o*c „ 



t < t c 
t c <t < t c + T c 
t c + T c <t 



We now assume that the control Liouvillian scales in- 
versely with the pulse duration C c — — with k c being a 
dimensionlcss super-operator of Lindblad form. Letting 
the pulse duration vanish r c — > , we approximate a <5- 
kick via £(t) = Cq + n c 5{t — t c ) (where Six) represents 
the Dirac-c) distribution), such that the solution reads 



Pit) 



e Cot p 

^Co{t — t c ) gK C gCot c 



t < t c 

tc < t 



(Bl) 



Appendix B: <5-kick propagators 

There are multiple derivations of delta-kick propaga- 
tors, we provide a simple pedestrians derivation based on 
piecewise constant time dependencies as used in scheme 
I. When £(i) is a matrix that does not commute with 
itself at different times, the equation p(t) = £(t)p(t) is 
impossible to solve analytically in the general case. How- 
ever, when the Liouvillian £(£) is piecewise constant, the 
solution may be readily obtained by conventionally prop- 
agating for a time period where the Liouvillian is con- 
stant, and using the resulting state as an initial value 
for the next propagation period with another constant 



that is, the constant propagator evolution is simply in- 
terrupted by the control operation e Kc at control time 

tc- 

Experimentally, it may be more reasonable to discuss 
smooth dependencies £(t) = Cq + C c it — t c ) with the 
assumptions C c ix < — t c /2) — C c ix > +t c /2) = and 
that J C c it)dt is independent of the pulse duration r c . 
The latter condition also implies that the maximum pulse 
height must scale inversely with the pulse duration. In 
the limit of vanishing pulse duration r c — > such a con- 
trol operation would also converge to a Dirac-J distribu- 
tion, such that during control, the baseline Liouvillian Cq 
may be neglected and we would have an effective propa- 
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FIG. 3: (Color Online) Comparison of a single (thin red curve 
with jumps, same realization in all panels) and the average 
of 100 (medium thickness, green) and 10000 (bold smooth 
curve, turquoise) trajectories with the solution from the ef- 
fective feedback master equation [Eq. (3), thin black] for the 
dot occupation (top), the number of particles on the left (mid- 
dle), and the number of particles on the right (bottom). The 
average of the trajectories converges to the effective feedback 
master equation result. The reference curve without feed- 
back (dashed orange) may be obtained from Eq. (2) or by 
using vanishing feedback parameters and demonstrates that 
the direction of the current may actually be reversed via 
sufficiently strong feedback. Parameters: Tl = Fr = T, 
h = 0.45, f R = 0.55, 5f = 5g = 1.0, 5f = Sf = -10.0, 
and TAt = 0.01. 



gator 



+t c /2 



n c = lim / C c (t')dt' 

r c ^0 J 



(B2) 



-Tc/2 



which must be dimensionless and inherits the Lindblad 
form from C c (t). 



Appendix C: Justification of scheme Ila 

When ultrashort (in the sense discussed in Ap- 
pendix B) control operations on the SET tunneling rates 
are only to be applied immediately after an electron 
jumps into or out of the dot, the derivation of an effective 
master equation is a bit more involved. We start from 
the inverse Fourier transform of the evolution generated 
by Eq. (2) 



p(^.»H)(t + At) 
J^ nR \At) 



C LJ^pC(XL,XR.)^-in-L-XL-in R -XR 

4tt 2 



where pit) denotes the unconditional dot density matrix 
at time t, and n^m the number of electrons that have 
tunneled to left/right reservoirs during the timestep At. 
For small At it suffices to consider single-particle jumps 
only. Expanding the propagator for small At and using 
the orthonormality relation f_ w e l ^ n ~ m ^ x dx = 27r<5 5 
obtain the conditional evolution equations 



we 



p<°'°>(i + At) = 

p(+i.°)(i + A*) = 

p^°\t + At) = 

p(°> +1 \t + At) = 

p^Ht + At) = 



[1 + At (T L Fl 
T L AtT+p(t) , 
T L AtTlp{t) , 
T R AtT+p(t) , 
T R AtT R p(t) , 



(C2) 



and the respective probabilities are given by the trace 
of these operators. To derive a master equation ac- 
counting for the average evolution of observables without 
feedback one may simply compute the weighted average 
p(t + At) = pM + + pt- 1 - ) + p( > +1 ) + p' -- 1 ), 

which would after solving for the finite difference scheme 
{pit + At) — p(t)) /At yield the original Liouvillian in 
Eq. (1) when At — > 0. In contrast, with (5-kick feedback 
operations (compare appendix B) e.g. 



+r/2 

Kj = J Cf ntm \t')dt' 

-T c /2 

= 5{F L {Q) + 5 R T R (Q) , 



(C3) 



compare Eq. (2) and similarly for no, where S a 

+t/2 

J Tj (t)dt' are dimensionless feedback parameters 

-T c /2 

(characterizing the product of height and width of the 

time-dependent control tunneling rates vH °{t)), we have 
to perform the average after applying the matching feed- 
back operation [in (I) or out (O)] 

p(t + At) = p M +e «o p (.+l,0) +e n lp (-l,0) 

+e «o p (0,+l) +e Ki p (0,-l) ( C4 ) 

After converting this into a finite-difference equation we 
would obtain the generator in Eq. (6) without count- 
ing fields. Now it is essential that during the control 
operations, electrons may tunnel through the junctions: 
Although the pulse duration is infmitcsimally small, the 
product of pulse height and pulse width remains con- 
stant, such that the tunneling probability is character- 
ized by the feedback strength. Therefore, counting fields 
may be re-introduced via a conditional master equation 
as was done when deriving Eq. (2) from Eq. (1). Finally, 
this yields the form of Eq. (6). 

For the trajectories we constrain ourselves to opening 
only one junction at a time for all control operations, e.g., 
8 L = 6 R = in order to avoid an unbounded number of 
(CI) particles tunneling during control. It is straightforward 
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to compute the trajectories as described in appendix A, 
sec Fig. 4. The computed trajectories do well converge to 
the results of the effective feedback master equation. The 
scheme is non- recursive, i.e., electrons tunneling during 
control operations do not trigger a further control oper- 
ation. 



FIG. 4: (Color Online) Comparison of a single and the av- 
erage of 100 and 10000 trajectories with the solution from 
the effective feedback master equation (6) for the dot occupa- 
tion (top), the number of particles on the left (middle), and 
the number of particles on the right (bottom). The average 
of the trajectories converges to the effective feedback master 
equation result. The single trajectory only records the state 
after all control operations have been performed, i.e., jumps 
with no net effect are not visible. Therefore, control opera- 
tions (bold kinks in top graph) are not always correlated with 
a net change in the occupation or the particle number left 
and right. The second last top kink for example stands for a 
out-control operation e K ° that was triggered by an electron 
jumping to the left (an information which is not provided by 
the QPC) contact. During the control operation, an electron 
jumped back from the left contact to the system, such that 
no net change in ridot, ni, nor tir occurred. Color coding 
and parameters are as in Fig. 3 where applicable. Feedback 
parameters were chosen as <5£ = <5§ = and S R — <5£ = 1.0. 



Appendix D: Justification of scheme lib 



In order to resum the control operations when applied 
recursively, it is essential that only one junction is opened 
at a time, formally expressed by choosing 5^ = S R = 
0. This implies that once a control operation has been 
triggered initially, the transport becomes unidirectional, 
which allows for a simple analytic rcsummation. 



J 



We note that in this case the control operations only depend on a single counting held and we have in their matrix 
exponential the simple decomposition 



ko(xl,Xr) 



V°( X l) =K + P8e +lXL + Vfe~ 1 ^ 



1 - (1 - )f L e+*» (1 - e" 5 - )(1 - f L ) 
B -ix*(l_ e -tf)/ L i_(i_ e -*2)(l-/ L ) 



Vie 



I„-iXR 



1 - (1 - e- 5 *)f R e+*«(l - e- A «)(l - f R ) 
e-««(l - e- s 'n)f R 1 - (1 - e-**)(l - f R ) 



(Dl) 



where the upper index labels the trigger process and the lower indices mark the parts responsible for no particle 
change (N), a particle jumping out of the system (O, possible for kj only to the right), and a particle jumping into 
the system (I, possible for kq only from the left). 



The fact that out of an empty system, no particle may 
jump out and vice versa for a filled SET system no further 
particle may jump in is formally reflected in the relations 
pOj-+ = v o j:+ = o and pijr- = V\T R = as well as 



L 

0. Therefore, the first-order feedback 
master equation for scheme Ha in Eq. (6) reduces in this 



case (Si = <jg = 0) to 



c 



(1) 

fb 



+ (P° 



(D2) 
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To generalize this scheme to recursive feedback scheme 
(lib) with a potentially infinite recursion depth requires 

to recursively follow the scheme £^ — > Cp^ 1 ^ by per- 
forming in each iteration the replacements 



VfA 

r' N A 



{Vlt + Vhe+^VfA, 
VnA, 

{V^ + Vfe-^V'oA 



(D3) 



for arbitrary operators A (in this recipe we have also used 
that VjVf = TgVjj = 0), which eventually leads to 



C 



— U^r,+1 R-Tr 



1'b 



lO 



+{T D N + 'P I N V?e-^) x 

oo 

V"V"Pq ■p0^n e +in(xR-XL) 



n=0 



+{V I N +V%V I e+^ R ) x 



+in{xR-Xh) 



x(T L e-^Fl +T R e- l *»Fx). (D4) 

After summing the von Neumann operator series we fi- 
nally obtain Eq. (7). 

Similarly, we may recursively construct the trajecto- 
ries numerically, see Fig. 5. The tunneling during mu- 
tually calling control operations is only stopped after no 
particle tunnels. At infinite bias and infinite feedback 
strength, the corresponding halting probability vanishes, 
which leads to a feedback catastrophe. Precursors of this 
are already observed at finite feedback strength and infi- 
nite bias in the large current and Fano factor - compare 
right columns ins tables I and II - or at infinite feed- 
back strength and finite bias in the current - compare 
the dotted curve in Fig. 2. 



Appendix E: Unconditional control leads to 
transport with the bias 

This statement is mainly based on detailed balance, 
where we assume without loss of generality /l < fn at 
the SET level. 

By evaluating the stationary states p = (1 — 
SET ) T of Liouvillians C, £e, and Cf one can 
see that their stationary occupation n T is within the 
transport corridor Jl < n SET < f R with the actual 
position depending on the respective tunneling rates. 
With any state within the transport corridor as ini- 
tial condition one can show that the propagators V 6 



n SET ,n' 



e CAt^ e c E At^ e c R At 



{ 

Pi+i 



(1 



,SET „SETaT 
i+1 ' ' 



, e Kl , e K ° } again map to occupations 
Ti T ) T = m nf ET ,nf ET f = V Pi 



within the corridor, i.e., for all fz < nf ET < fa we have 
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FIG. 5: (Color Online) Comparison of a single (only net 
changes as in Fig. 4) and the average of 100 and 10000 tra- 
jectories with the solution from the effective feedback master 
equation (7) for the dot occupation (top), the number of par- 
ticles on the left (middle), and the number of particles on 
the right (bottom). The average of the trajectories converges 
to the effective feedback master equation result. Electrons 
jumping during control operations may now recursively trig- 
ger further control operations - marked by the height of the 
bold kinks in the top graph. Color coding and parameters are 
as in Fig. 4. 



/l < n ifi < Ir- This implies that once any uncon- 
ditional iterative control scheme has entered the trans- 
port corridor (large times), it will not be able to leave 
it again, regardless of the timesteps At, control parame- 
ters, baseline tunneling rates, and the order of the opera- 
tions. In contrast, with the conditioned feedback scheme 
one is always outside this corridor with n^ ET s {0, 1}. 
With inserting counting fields at e.g., the right junction 
V — > V(xr) one may now calculate the mean particle 
number tunneling during an iteration step 

An t = (-i)d XR Tr{r( XR )(l~n^,n^f}\JEl) 



when the initial state is within the transport corridor 
(/l < n f ET < Jr.)- The outcome is that in this case the 
electron current always points from right to left when 
/l < Ir, i-C-, for similar temperatures with the bias. 

Note also that for unconditional switching between the 
Liouvillians (scheme I), this is not a Parrondo game any- 
more, since a third possibility of game outcome (no tun- 
neling event at all) is included. 
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Appendix F: Nonlinear Current- Voltage Characteristics 



The stationary current (e.g., at the right junction) may either be calculated from the cumulant-generating function 
or via the relation / = (— i)Tr |<9 XB C(xl,Xr)\q q p|, where £(0, 0)p = and reads for scheme I 



iL = 



r L e d £ (1 - f L ) + T R e 6 n(l - f R ) + T L e"tf L + T R e d nf f 



(Fl) 



where insertion of the Fermi functions f a = [e^°( c_A1 °) + 1] _1 at similar temperatures /3l = fl R = f3 and symmetric 
chemical potentials /i£ = +V/2 and fi R = —V/2 yields the full nonlinear current-voltage characteristics displayed in 
Fig. 2. Naturally, the zero-feedback case 6f = 8^ = 5^ = 8^ = reproduces the known results. 

For feedback schemes Ha and lib the expressions become a bit lengthy, such that we only give the maximum 
feedback limit 8^ = S R = and S R = 82 — > oo for the current 



jlla 



- f R ) [r|(i - f L ) 2 + ryi] + t l t r [f L -fl-f 2 R - SlSh - Vlfl + VlUIl + f R j\ 
r L [i - /l(i - f L + 1 - f R )} + r R [l - (l - f R )(f L + f R )] 



T m = r|/L(i - h)(i - f R ) + ry L / fl (i - f R ) + T L T R [f L - fl - fjf R + f L p R ] 

lh ' T L [l-f L (l-f L + l-f R )]+r R [l-(l-f R )(J L + f R )] ■ ( ' 

Again, insertion of the Fermi functions at similar temperatures and symmetric chemical potentials yields the nonlinear 
current-voltage characteristics in Fig. 2. It also becomes evident that at reverse infinite bias (/^ — > and f R — > 1), 
feedback with 8 ! L = S R = cannot overcome the large bias voltage, such that the same result as without feedback is 
obtained. Inserting these highly nonlinear dependencies in the differential equation V = eI(V)/C finally yields the 
inset of Fig. 2. 



